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Using a transport model coupled with a phase-space coalescence after-burner we study the triton- 
''He ratio, relative and differential transverse flows in semi-central ^^^Sn -i-^^* Sn reactions at a 
beam energy of 400 MeV/nucleon. The neutron-proton ratio, relative and differential flows are 
also discussed as a reference. We find that similar to the neutron-proton pairs the triton-'^He pairs 
also carry interesting information about the density dependence of the nuclear symmetry energy. 
Moreover, the nuclear symmetry energy affects more strongly the t-^He relative and differential flows 
than the -k~ /tv^ ratio in the same reaction. The t-^He relative flow can be used as a particularly 
powerful probe of the high-density behavior of the nuclear symmetry energy. 

PACS numbers: 25.70.-z, 25.70.Pq, 21.65.Ef 
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I. INTRODUCTION 

The density dependence of nuclear symmetry en- 
ergy especially at supra-saturation densities is among 
the most uncertain properties of neutron-rich nuclear 
matter However, it is very important for nuclear 
structureUll], heavy- ion reactions 0, d, 0, H, H, and 
man y p henomena/processes in astrophysics and cosmol- 
ogy [111 . fT3 . [ist . Heavy-ion reactions especially those in- 
duced by radioactive beams provide a unique opportunity 
to constrain the symmetry energy at supra-saturation 
densities in terrestrial laboratories. Various probes using 
heavy-ion reactions have been proposed in the literature, 
see, e.g., ref.[l^ for the most recent review. It is partic- 
ularly interesting to mention that, besides many signifi- 
cant results about the symmetry energy at sub-saturation 
densities, see, e.g., refs.B 111 M, 0, H E, US, 
circumstantial evidence for a rather soft symmetry en- 
ergy at supra-saturation densities has been reported very 
recently[21| based on the IBUU04 transport model|23| 
analysis of the 7r~/7r+ data taken by the FOPI Collabo- 
ration at SIS/GSl[23|. However, many interesting issues 
remain to be resolved. Thus, to constrain tightly and re- 
liably the nuclear symmetry energy especially at supra- 
saturation densities, much more efforts by both the nu- 
clear physics and the astrophysics communities are still 
needed. 

In the present work, we have the following two main 
purposes. Firstly, there is an urgent need to verify the 
conclusion about the soft symmetry energy at supra- 
saturation densities required to reproduce the FOPI 
■K~ j-K^ data within transport model analyses pli [2^. 
It is better if this test can be done with not only 
more 'k~ jix^ data but also other sensitive obscrvablcs 



in the most neutron-rich reactions possible. We will thus 
make predictions using the same transport mo del for 
doing this test. Secondly, it was predicted that the 
neutron-proton differential flow is another sensitive probe 
of the hig h-density behavior of the nuclear symmetry 
energy |24|. However, it is difficult to measure observ- 
ables involving neutrons. One question often asked by 
some experimentalists is whether the triton-'^He pair may 
carry the same information as the neutron-proton one. 
We will try to answer this question quantitatively by 
coupling the IBUU04 calculations to a phase-space co- 
alescence after-burner. Indeed, we found that, similar to 
the neutron-proton pair, the triton-'^He relative and dif- 
ferential transverse flows are sensitive to the high-density 
behavior of the nuclear symmetry energy. They can be 
used to test indications about the high-density behavior 
of the symmetry energy observed earlier from analyzing 
the 7r~/7r+ data. 



II. SUMMARY OF THEORETICAL MODELS 

Our study is carried out based on the IBUU04 version 
of an isospin and momentum dependent transport model 
and the simplest phase-space coalescence after-burner. 
For completeness and consistency we recall here a few 
major features of the IBUU04 transport model most rel- 
evant to the present study. More details of the model can 
be found in Refs. [13] ■ The single nucleon potential is 
one of the most important inputs to BUU-like transport 
models for nuclear reactions. In the IBUU04 transport 
model, we use a single nucleon potential derived within 
the Hartree-Fock approach using a modified Gogny effec- 
tive interaction (MDI) [2^, i.e.. 
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Here (5 = (p„ — pp)/p is the isospin asymmetry of the nu- 
clear medium. In the above r = 1/2 (—1/2) for neutrons 
(protons) and t ^ t' ; a — 4/3; fri'r'iP) is the phase space 
distribution function at coordinate f and momentum p. 
The parameter x was introduced to mimic predictions on 
the density dependence of symmetry energy Egymip) by 
microscopic and/or phenomenological many-body theo- 
ries. The parameters Au{x) and Ai(x) depend on the x 
parameter according to 

2_S 2_B 

Au{x) = -95.98 -a; -, Ai{x) = -120.57 + x -. 

cr + 1 (T + 1 

(2) 

The coefficients in Au{x) and Ai{x) and the parame- 
ters B ,Cr,T,Cr,T' and A were obtained by fitting the 
momentum-dependence of the U {p, S, p, r, x) to that pre- 
dicted by the Gogny Hartree-Fock and/or the Brueckner- 
Hartrcc-Fock (BHF) calculations, the saturation proper- 
ties of symmetric nuclear matter and the symmetry en- 
ergy of about 31.6 MeV at normal nuclear matter density 
Po ~ 0.16 fm~^. The incomprcssibility Kq of symmetric 
nuclear matter at po is set to be 211 MeV consistent with 
the latest conclusion from studying giant resonances [2^ . 



The last two terms in Eq. ([T]) contain the momentum- 
dependence of the single-particle potential. The momen- 
tum dependence of the symmetry potential stems from 
the different interaction strength parameters Ct,t' and 
Ct,t for a nucleon of isospin r interacting, respectively, 
with unlike and like nucleons in the background fields. 
More specifically, we use Cu = CunUke = —103.4 MeV 
and Ci = Ciike = —11.7 MeV. With these parame- 
ters, the isoscalar potential estimated from {Uneutron + 
Uproton)/'^ agrees reasonably well with predictions from 
the variational many-body theory [l^, [1^ • At the nor- 
mal nuclear matter density po, it is consistent with 
the isoscalar nucleon optical obtained from the nucleon- 
nucleus scattering experiments [2^. The strengthes of 
the corresponding isovector (symmetry) potential can be 
estimated from (C7„ — Up)/26. At po, by design, the sym- 
metry potential is independent of x and is consistent with 
the Lane potential extracted from the nuclcon-nucleus 
scattering experiments and the (p,n) charge exchange re- 
actions pi| . 

The corresponding MDI symmetry energy can be writ- 
ten as [13 
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where Pf = h{i-K^ ^Y^^ is the Fermi momentum for sym- 
metric nuclear matter at density p. 

The IBUU04 model can use either the experimental 
nucleon-nucleon (NN) cross sections or the in-medium 
NN cross sections calculated using an effective-mass scal- 
ing model consistent with the single particle potential 
used |15j . In the present work the in- medium NN cross 



sections are used. The isospin-depcndent Pauli blocking 
has been implemented by evolving and checking explic- 
itly neutron and proton phase-space distributions sep- 
arately. The coordinates of nucleons in the colliding 
nuclei are initialized randomly according to the neu- 
tron/proton density profiles predicted by the Skyrme- 
Hartree-Fock approach. The corresponding Fermi mo- 
menta are calculated using the local Thomas-Fermi ap- 



proximation. The initial state generated in such way 
is rather stable for several hundred fm/c without ap- 
preciable particle emission in evolving a single nucleus 
with a momentum-independent mean-field. However, as 
it is widely known, see, e.g., ref, [2^, momentum de- 
pendent mean-field makes the initial state less stable. 
For instance, in evolving a single ^^^s^ nucleus with the 
IBUU04, at 40 fm/c (by which most of the nucleons 
should have freezed-out in heavy-ion reactions at a beam 
energy of 400 MeV/nucleon) the average ratios of emit- 
ted/initial protons and neutrons are 0.6/50 and 1.7/74, 
respectively, in calculations using 600 test-particles per 
nucleon. The ratios go up to about 1.5/50 and 3.4/74, 
respectively, in calculations using 200 test-particles per 
nucleon. Our results presented in the following are ob- 
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tained using totally 10,000 events in each case by using 
200 test-particles per nucleon in each run of the simula- 
tion. 

Because most BUU-type transport models including 
the IBUU04 are incapable of forming dynamically realis- 
tic nuclear fragments, some types of after-burners, such 
as statistical and coalescence models, are normally used 
as a remedy. This kind of hybrid models can be used 
to study reasonably well, for instance, nuclear multifrag- 
mentation, see, e.g., ref. [31I. [331 , collective flow o: ' 

light fragments, see e.g., SsT 34|, l35j and the formation 
of hvpernuclei [36l| . There are, however, some remain- 
ing issues, such as the freeze-out time of fragments that 
is related to the time of coupling the transport model 
with the after-burner, etc. There are also interesting 
work in using advanced coalescence models[33, [111, see, 
e.g., refs.[39l.l40j. We notice here that, several advanced 
cluster recognition routines, such as, the Early Cluster 
Recognition Algorithm (ECRA) , the Simulated An- 
nealing Clusterization Algorithm (SAC A) [S^l, have been 
put forward in recent years. For the purposes of the 
present exploration, however, we use the sirnplest phase- 
space coalescence model, see, e.g., refs.fs^. [35|. where a 
physical fragment is formed as a cluster of nucleons with 
relative momenta smaller than Pq and relative distances 
smaller than Rq. The results presented in the follow- 
ing are obtained with Pq = 263 MeV/c and i?o = 3 
fm. This simple choice may thus limit the scope and 
importance of our study here. For instance, we shall 
limit ourselves to studies of the relative/differential ob- 
servables for neutron-proton and t-'^He pairs without at- 
tempting to study pairs of the heavier mirror nuclei. An 
extended study including the heavy mirror nuclei using 
the advanced coalescence and/or earlier cluster recogni- 
tion methods is planned. 



III. RESULTS AND DISCUSSIONS 

Noticing that the FOPI 7r^/7r+ data favors the sym- 
metry energy with x = 1 at supra-saturation densities [21 1 
while the NSCL/MSU isospin diffusion data favors a sym- 
metry energy at sub-saturation densities between those 
with X = and x = — l[ll, [I!], for comparisons we use 
here a; = 1 and 2; = — 1 as two limits of the symmetry 
energy at high densities. The corresponding symmetry 
energy functionals are depicted in the inset of Fig. [T] 
They represent a typically stiff {x = —1) and a very soft 
{x = 1) symmetry energy at supra-saturation densities. 

In the FOPI experiments [23j , the 7r~/7r+ ratio was 
measured down to a beam energy of 400 MeV/nucleon 
where it shows the largest sensitivity to the high-density 
behavior of the nuclear symmetry energy pl|. The maxi- 
mum density reached in central Au-|-Au reactions at 400 
MeV/nucleon is about 2.5po. It is well known that pions 
are most abundantly produced in the central collisions. 
However, the 7r~/7r"'" ratio is almost a constant from most 
central to mid-central impact parameters [2 1|. It is also 
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FIG. 1: Multiplicities of triton and ''He as a function of the 
reduced CM. rapidity in the reaction of ^'^'^ Sn -\-^'^'^ Sn at a 
beam energy of 400 MeV/nucleon and an impact parameter 
of 5 fm with the symmetry energy parameter x = 1 (soft) and 
X = —1 (stiff), respectively. 



well known that transverse collective flow is zero in head- 
on and grazing collisions but is the largest in mid-central 
collisions. Considering all of the above, as an example, 
we study ^^^Sn -|-i24 reaction at an incident beam en- 
ergy of 400 MeV/nucleon and an impact parameter of 5 
fm. This reaction will be available at FAIR/GSI in the 
near future. We note that the maximum density reached 
in this reaction is about 2po- 

Before we study the triton-'^He relative and differential 
flows, it is necessary to first examine the yields of tri- 
ton and '^He and their ratio. Although these observables 
and their dependence on the symmetry energy have been 
studied before especially at lower incident energies, see, 
e.g., ref. m, HO], they serve as useful references for mea- 
suring the symmetry energy effects on the 7r~/7r"'" ratio 
and the flow observables. Shown in Fig.[T]are the rapidity 
distributions of triton and '^He for the ^^^Sn -(-124 j.g_ 
action at an incident energy of 400 MeV/nucleon and an 
impact parameter of 5 fm. It is seen that these clusters 
are produced mostly at mid-rapidities from the partici- 
pant region. The yields, however, are not sensitive to the 
symmetry energy. This is not surprising. Even at much 
lower energies where effects of the symmetry energy is 
stronger, the neutron and proton yields themselves are 
not so sensitive to the symmetry energy as the yields are 
dominated by the isoscalar part of the nuclear mean field 
and nucleon- nucleon collisions. 

At Fermi energies, the ratio of neutrons to protons 
or that of mirror nuclei were shown to carry more 
information about the symmetry energy as effects of 
the isoscalar potential can be largely cancelled in the 
ratiosjH [H, SO, [H, HI]. However, at significantly 
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FIG. 2: Comparisons of particle ratios in the same reaction 
as in Fig. [T] The horizonal line at 1.56 (2.43) in the left 
(right) window is the reaction system's neutron/proton ratio 
(squared) . 



higher energies, e.g., 400 McV/nuclcon, except for high 
transverse momentum particles especially those being 
squeezed-out perpendicular to the reaction plane [45j. the 
neutron to proton ratio becomes less sensitive to the 
symmetry energy than at lower beam energies [23, l46j . 
For comparisons, the t/^Hc ratio together with the free 
and bound neutron to proton rations are shown in the 
left window of Fig. [5] for the ^^-^Sn -|-^^^ Sn reaction. 
Here, the bound neutrons and protons are the nucle- 
ons in the fragments with A > 2. The horizonal line 
at 1.56 is the reaction system's neutron/proton ratio 
(Nt + Np)/(Zt + Zp). Several interesting observations 
can be made here. Firstly, the (n/p)bound (("■/p)frcc) is 
significantly less (higher) than the neutron/proton ratio 
of the reaction system. This is the well known isospin 
fractionation phenomenon [i^, IH, [4^ which is reduced 
here by the production of charged pions. Moreover, at 
mid-rapidity the (n/p) bound shows appreciable sensitivity 
to the variation of the symmetry energy. Unfortunately, 
both the {n/p)iree and t/^He ratios show very little sen- 
sitivity to the variation of the symmetry energy within 
statistical error bars, except around the projectile and 
target rapidities of y/j/beam = ±0.5. 

It is especially worth noting that the t/'^He ratio is 
much higher than the free and bound neutron/proton 
ratios. A few more comments about this observa- 
tion are in order here. First of all, we notice that 
the assumption of t/^He={n/p)ircc and that they are 
equal to the (Nx 4- Np)/(Zt + Zp) have been widely 
used in momentum-space coalescence models in the lit- 
erature, especially at high energies. Our results here 
and also those reported earlier from transport model 
simulations [U |50| | indicate that the coordinate-space 
correlation is important in considering the cluster for- 
mation. This feature has also been observed in the par- 
ton hadronization process in ultra-rclativistic heavy-ion 



collisions by the quark coalescence model (5l|. In fact, 
a similarly high t/^He ratio has also been observed in 
the experimental data of heavy-ion reactions at Fermi 
energies, see, e.g., ref.[3l|. The possible explanations for 
the high n/p and t/^He ratios include the isospin frac- 
tionation and the overwhelmingly preferential production 
of symmetric light clusters such as deuteron and alpha 
particlesjs^l in these reactions. Our analyses by turn- 
ing on/off the coalescence after-burner, the Coulomb and 
symmetry potentials indicate that both mechanism are at 
work. It is also interesting to mentioning that, at Fermi 
energies, using a freeze-out temperature extracted from 
the experiments the rather high t/'^He ratio can also be 
well reproduced within a statistical model [l^l assuming 
that the difference in the chemical potentials of triton 
and ''He is dominated by their Coulomb potentials |31i]. 

The rapidity distributions of the tt^/tt^ ratio and the 
bound protons are shown in the right window of Fig. [21 
Interestingly, comparing the ratios of all particle-pairs 
shown in both windows of Fig. [51 it is obvious that the 
7r~/7r+ ratio shows the highest sensitivity to the symme- 
try energy. More quantitatively, effects of the symmetry 
energy on the ratio of the total yields of charged pions 
is about 10% by varying the x parameter from —1 to 
1. This is less than the approximately 20% effect ob- 
served in the head-on collisions between two Au nuclei 
at the same beam energy. This is probably because of 
the significantly smaller size in ^^^Sn +124 although 



it is more neutron- rich 5J| . Moreover, by comparing the 
rapidity distributions of the 7r~ /tt^ ratio and the bound 
protons one sees clearly the well-known Coulomb focus- 
ing effects on the 7r~/7r+ ratio, namely, more 7r~'s (tt+'s) 
are attracted (repelled) towards (away) from the target 
and projectile residues [ssllSGj where most of the protons 
are located in the semi-central reactions considered. 

We now investigate whether the transverse collective 
flows of triton and "'He can be used to probe the symme- 
try energy. Firstly, we examine in Fig. [31 their transverse 
flows individually. The average CM. transverse momen- 
tum per nucleon < Px/A > in the reaction plane is de- 
fined as 



<PxlA> (y) 



1 



N{y) 



Niy) 



(4) 



where N{y) is the total number of fragments of mass A 
in the rapidity bin at y. The correlation between the 
< PxJA > and rapidity y reveals the transverse collective 
flow [531 ■ It is seen that '^He clusters show a stronger flow 
than triton clusters. This is mainly due to the stronger 
Coulomb force experienced by the '^He clusters. More 
interestingly, the transverse flow of ^He clusters show 
appreciable sensitivity to the variation of the symmetry 
energy. 

The transverse flow is a result of actions of several fac- 
tors including the isoscalar, symmetry and Coulomb po- 
tentials and nucleon-nucleon scatterings. It is well known 
that the transverse flow is sensitive to the isoscalar poten- 
tial. Given the remaining uncertainties associated with 
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FIG. 3: Triton and ^He transverse flows (in unit of MeV) as 
functions of the reduced CM. rapidity in the same reaction 
as in Fig. [T] 



the isoscalar potential and the small size of the symmetry 
energy effects, it would be very difficult to extract any re- 
liable information about the symmetry energy from the 
individual flows of triton and '^He clusters. Thus tech- 
niques of reducing effects of the isoscalar potential while 
enhancing effects of the isovector potential are very help- 
ful Hi, SI, m, m, [H. We thus study in Fig. [4] the 
triton-'^He relative and differential flows. The relative 
flow is given as 




1 ^* 

< v'JA >-< pI'^VA >= j^Y^^JA 



The triton-^He differential flow reads 
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FIG. 4: Triton-'^He relative and differential flows (in unit of 
MeV) as a function of the reduced C.M. rapidity in the same 
reaction as in Fie. [T] 
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where Nt, Na-^^ are the number of triton and "^He in the 
rapidity bin at y. From the upper panel of Fig. [H it is 
seen that the triton-'^He relative flow is very sensitive to 
the symmetry energy. Because of the larger slope of the 
^He flow, the triton-'^He relative flow shows a negative 
slope at mid-rapidity. Effects of the symmetry energy on 
the differential flow shown in the lower panel, however, is 
relatively small. Although the "^He flow is more sensitive 
to the symmetry energy, the small number of "^He clus- 
ters (as shown in Fig. [T|) makes the '^He flow contributes 
less to the triton-'^He differential flow (as indicated in Eq. 
(O). The triton-'^He differential flow is therefore domi- 
nated by triton clusters. Consequently, it is less sensi- 
tive to the symmetry energy than the triton-'^He relative 



flow. The slope F{x) = d < Px/A > /d{y/ybearn) of the 
transverse flow at mid-rapidity can be used to character- 
ize more quantitative the symmetry energy effects. We 
found that for the t-'^Hc relative flow, F{x = I) « —74 
MeV/c and F{x = -I) « -22 McV/c, respectively. For 
the t-^Hc differential flow, F{x = I) « 21 MeV/c and 
F{x ~ —1) w 42 MeV/c, respectively. Compared to 
the tt^/tt^ ratio in the same reaction, the symmetry en- 
ergy effects on the t-'^He relative and differential flows are 
much stronger. Thus, especially the t-'^He relative flow 
can be used as a very useful and independent tool to test 
the soft symmetry energy at supra-saturation densities 
extracted from studying the 7r^/7r+ ratio 21 1. 



For comparisons, we now study the relative and differ- 
ential flows for free neutron-proton pairs in Fig. [S] ft is 
seen that they have the same features as the relative and 
differential flows for triton-"^ He pairs. The larger sym- 
metry energy effects at positive rapidities are due to the 
more neutron-rich of projectile. More quantitatively, in 
terms of the slope parameter F, for the n-p relative flow, 
F{x = 1) w -53 MeV/c and F{x ^ -I) w -25 MeV/c, 
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FIG. 5: Neutron-proton relative and differential flows (in unit 
of MeV) as a function of the reduced CM. rapidity in the same 
reaction as in Fig. [T] 



respectively. For the n-p differential flow, F{x = 1) w 20 
MeV/c and F{x = —1) w 36 MeV/c, respectively. Com- 



paring the results in Figs. [5] and [4] and the slope pa- 
rameters F for t-'^He and neutron-proton pairs, we can 
conclude that they are almost equally useful for probing 
the density dependence of the nuclear symmetry energy. 



IV. SUMMARY 

In summary, using a hybrid approach coupling the 
transport model IBUU04 to a phase-space coalescence 
after-burner we studied the t-^He relative and differen- 
tial flows in semi-central ^^"^Sn -|-124 reactions at an 
incident energy of 400 McV/nucleon. We found that the 
nuclear symmetry energy affects more strongly the t-'^He 
relative and differential flows than the tt^/tt^ ratio in 
the same reaction. The t-'^He relative flow can be used 
as a particular powerful probe of the high-density be- 
havior of the nuclear symmetry energy. It can be used 
to test the indications about the symmetry energy at 
supra-saturation densities observed in the analysis of the 
■K~ /t:^ data from heavy- ion reactions. 
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